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Abstract 

We calculate the two-loop matching corrections for the gauge couplings at the 
Grand Unification scale in a general framework that aims at making as few assump- 
tions on the underlying Grand Unified Theory (GUT) as possible. In this paper we 
present an intermediate result that is general enough to be applied to the Georgi- 
Glashow SU (5) as a "toy model" . The numerical effects in this theory are found to 
be larger than the current experimental uncertainty on a s . Furthermore, we give 
many technical details regarding renormalization procedure, tadpole terms, gauge 
fixing and the treatment of group theory factors, which is useful preparative work 
for the extension of the calculation to supersymmetric GUTs. 

PACS numbers: 12.10.Kt 12.38.Bx 02.20.Sv 11.25.Db 



1 Introduction 



Grand Unified Theories (GUTs) provide an appealing framework for physics beyond the 
Standard Model (SM) of particle physics. Particularly supersymmetric (SUSY) GUTs 
have gained a lot of interest in the past decades as they seem to be consistent with the 
measured values of a s (Mz), a em {Mz) and sin9w(Mz) [IH3] and offer many other beau- 
tiful features. Up to date most Renormalization Group (RG) analyses that study gauge 
coupling unification use two-loop Renormalization Group Equations (RGEs) and one-loop 
matching at the SUSY and GUT scales [H4TT] . These matching corrections arise from in- 
tegrating out heavy particles at the SUSY and GUT thresholds. They depend sensitively 
on the mass splittings between the heavy particles and can be used to constrain the mass 
spectrum of the theory. As experimental accuracy is increasing, also higher order correc- 
tions become more and more important. It has been shown that O(al) effects at the SUSY 
decoupling scale can be as large as the current experimental uncertainties on the gauge 
couplings [T2TTH] . Furthermore, in GUT models that contain large representations, as for 
example the so-called Missing Doublet Model [T5|[TB]. the decoupling scale dependence at 
the GUT scale can exceed the experimental uncertainty by an order of magnitude [T%] . 

These facts encourage us to try to improve the theoretical accuracy of the unification 
study and aim at a complete three-loop RGE analysis. This requires three-loop RGEs 
for the SM, the Minimal Supersymmetric Standard Model (MSSM) and the SUSY GUT 
model under consideration and two-loop matching corrections at the SUSY scale and the 
GUT scale. In Quantum Chromodynamics (QCD) the gauge /3 function in the modified 
minimal substraction scheme (MS) is known to four loops [HI US], though a complete 
three-loop SM gauge function is still missing. For the MSSM full three-loop RGEs 
are available [191 120] in the so-called DR scheme. The same is true for the most general 
single gauge coupling theory in MS and a general supersymmetric GUT in DR [21TI22] . 
Matching corrections at the SUSY scale are known to one- loop order for the £7(1) and 
577(2) gauge couplings a±/2 [7J and to two-loop for the strong coupling a s [T2KI3]- For 
the GUT scale thresholds a general formula at the one-loop level is known [23H2HE3], 
but two-loop corrections have still been missing up to date and thus are not included in 
present unification analyses [H-[6l[T0 | [TT | [n] . 

The aim of this paper is to provide a first step towards a general formula for the two-loop 
GUT matching corrections in a similar fashion as the one-loop corrections in refs. [231 
[2HE3]. Unfortunately it turns out that at the two-loop level it is much harder to carry 
out the calculation in a general way, making as few assumptions about the underlying 
GUT model as possible. Therefore, the result given in this paper is not yet applicable to 
SUSY GUT models, as it is not yet general enough. Nevertheless, by applying it to the 
Georgi-Glashow 577(5) [M] as a "toy model", we provide an important intermediate step 
on the way to a full three-loop unification analysis. The generalization for SUSY GUTs 
actually is in progress. 
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The remainder of the paper is organized as follows: in the next section we describe the 
theoretical framework that is used for the calculation. There the Lagrangian is defined 
and it is shown how delicate issues, as gauge fixing and renormalization are carried out. 
In section [3] we present an academic study of our general formula when applied to the 
Georgi-Glashow model, which is the simplest (yet already ruled out) possible GUT model. 
In particular the issue of reducing the decoupling scale dependence is discussed. Finally, 
we present our conclusions. In appendix [A] of this paper we describe the procedure of 
defining and reducing the group theory factors that appear in our calculation. The rest 
of the appendix is dedicated to some supplementary material to the main text. 



2 Theoretical framework 

Since there is a host of well motivated GUT models, we want our final result to be applica- 
ble to as many of them as possible. Therefore, it is desirable to carry out the calculation 
of two-loop matching corrections at the GUT scale in a framework that makes as few 
assumptions on the underlying GUT model as possible. The idea is to have a general for- 
mula that depends on the Casimir invariants and the spectrum of the theory. Choosing a 
specific model specifies those Casimir invariants and gives an expression that depends only 
on the masses and couplings of the model. The actual calculation that is presented in this 
paper is not yet done in full generality, but makes some additional assumptions about the 
model. These assumptions are described in subsection 12.31 Nevertheless, we present the 
theoretical framework that is needed for the calculation of the relevant Green's functions 
(almost) as general as possible below in order to be armed for future improvements of the 
calculation. 

2.1 The Lagrangian 

We consider a general renormalizable quantum field theory defined by the following La- 
grangian: 

£ = --F» V F« V + tfiiyftf + l{p*$fD,& - U($) + C Y + £ gf + £ gh . (1) 

The chiral Dirac fermion field $ and the real scalar field $ reside in (not necessarily 
irreducible) representations of the gauge group G. The dynamics of the gauge field 
that transforms according to the adjoint representation of G, is described by the Yang 
Mills curl F£ u = d»A u a - d u A» + gf aM A^A v r Moreover, V($), Cy, £ gf and £ gh are the 
scalar potential, the Yukawa interactions, the gauge fixing and the ghost parts of the 
Lagrangian, respectively. They will be described in detail later in this section. V(<&) is 
chosen such that the scalar field $ = v + <&' contains one G-irreducible subspace that 
develops a vacuum expectation value (vev) v, that breaks G down to the SM gauge group 
H k G k = SU(3) x SU(2) x U(l). Models that have more than one vev of 0(M G ), where 
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Mq is the Grand Unification scale, are not covered by our framework yet. The indices 
a,/3, 7... belong to the adjoint representation and label the generators of G which fulfill 
the commutation relations 

[T a , T 13 } = if a ^ T 7 and also: [f a , f 13 } = if a ^ T 7 , (2) 

with the structure constants . We use the tilde to denote the generators of the 
reafl scalar representation which fulfill (T a ) T = —T a . The generator that acts on the 
fermion field satisfies (T a y = T a . Again, T a and T a need not necessarily be defined 
on irreducible representations of G, but can also have block diagonal form. In order to 
distinguish between broken and unbroken generators, we introduce the notation: 

{a} = + = {A} + {a} , (3) 

i i 

where Ai label the broken generators of G belonging to the SM-irreducible subspace 
labeled by i. If there is only one SM-irreducible subspace in the adjoint representation of 
G, as e.g. is practically the case in 577(5), we can omit the sub-index i. In contrast, 
label the unbroken generators belonging to the subgrou 

f^v = 0, 

t M v ^ 0. (4) 

The Lagrangian in eq. ([T|) is invariant under local gauge transformations with the real 
parameter 6 = 9{x): 

* -> m-ie a T a m, 

Al -+ A^ + r^Al--d,9 a . (5) 

The covariant derivatives are defined as: 

Dfl = (d,-igT a A«)y, 

= (d fl -igf a A$$. (6) 

Using eq. (|2J), eq. (BJ, and f aib i Ak = (cf. appendix |A]), the gauge-kinetic term for the 
scalar field $ = ?; + $' can be written as: 

1 -{D^) T D^ = ^{d^'fd^' + l -g 2 vf A f B vA» A A, B 

+igvf%$A ll A + ig 2 f ABa A B A iia vf A <$>' 

+ig®f a d^'A^ + ±g 2 &f a f f> &A£A, ia 

+g 2 vf A f B &A A \A^ B , (7) 



This is no loss of generality since every complex scalar can be written as two real scalars. 
2 Please note the different meanings of the sub-index i when attached to the capital adjoint index 
opposed to when attached to a lowercase adjoint index. 
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where we can identify the (diagonal) gauge boson mass matrix 



(M; 



X)AiBi 



g 2 vf Ai f B H 



with eigenvalues denoted by Mx r Again, the sub-index % labels the SM-irreducible sub- 
space that is meant, because each SM-irreducible subspace can be assigned to a definite 
gauge boson mass. Note that the position of the adjoint indices a,Ai,a{... is irrelevant. 
Furthermore, it is understood that a partial derivative acts only on the single field, which 
is next to it. The gauge- kinetic term for the scalars contains the undesired quadratic 
mixing igvT A d^A A between Goldstone bosons and heavy gauge bosons. As we will see 
in a moment, the gauge fixing Lagrangian £ g f can be chosen in such a way that this term 
is cancelled, at least at tree-level. 

In order to fix the gauge, we choose the gauge fixing functional [35] 



fAt 



d^ A -ig^i 2l vT A ^' 



(9) 



Note that we have chosen two distinct gauge parameters £1 and £2 for each SM-irreducible 
subspace of the heavy gauge bosons. They renormalize differently and thus can only be 
equated with each other after renormalization. Otherwise not all the Green's functions can 
be made finite. In the same way each SM group factor receives its own gauge parameter rji. 
This subtlety arises first at the two-loop level and is not relevant for computing one-loop 
matching coefficients. Although there are other ways to treat the gauge fixing [23 , 36 -38J , 
we find this one most convenient for our purposes. Employing the procedure described 
e.g. in ref. [39], eq. (|9]) gives us the gauge fixing Lagrangian and the ghost interactions: 



iVf 2 - - V f 2 



2£ 



— {d^ A f - ^g 2 £ 2i &f A *vvf A ^' + tgJ^vf A «S>'d»A A > 



li 



(10) 
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<h = Y,[ d " c AS 5 ^B 3 d,-gf AiBja A«)c Bj 

ij 

-g 2 Vb&i vT Ai T Bj v c A c Bj - g 2 y^^ vf A *f B ><S>'c j Ai c Bj 

~ [9(Vj/^ii) 1/4 fA i b j B i d fl c t A .A Bi c bj + ig 2 (£tv 1 / 4 £l{ 2 )fA i b j B i vf B ^'c Ai c b : 

ij 
ij 

+ E d " c U 5 ^A - 9fa^)c bi . (11) 

i 

Here and c ai denote the ghost fields belonging to the heavy and light gauge bosons, 
respectively. Note again that for an SU(5) GUT we could do the replacement A, — >■ A 
for the capital adjoint indices and the notation would become less clumsy. Here, however, 
we keep the sub-index i in order to stay as general as possible. As mentioned before, 
after partial integration the term ig^^i/iu vT^&d^A^ in eq. fTTOT) exactly cancels the 
corresponding term in eq. (J7|) at tree level, where £ij = ^ is a valid choice. However, when 
considering higher orders in perturbation theory, the bare gauge parameters £u and 
are not equal to each other and the above term must be kept explicitly as a counterterm in 
our calculation (cf. also subsection 12. 21) . The quadratic term in eq. ( flOl) can be identified 
with the (unphysical) Goldstone boson mass matrix: 

^coid = \g 2 Y,&i fAi ™ fAi ( 12 ) 

i 

with the property Tr(MQ old ) = Y2i M^. Df , where Df is the dimension of the i-ih 
SM-irreducible representation of the heavy gauge bosons. From the Goldstone theorem it 
follows that the structure v T Ai projects on the subspace of Goldstone bosons, i.e. on the 
subspace that obtains no mass term from V($). Hence, the matix Mq oM has non-zero 
entries only on the subspace of Goldstone bosons. 

For V($) we consider the most general renormalizable scalar potential with the discrete 
symmetry 

V($) = ~t4j^j + I k jk i $i®j®k®i , (13) 
with totally symmetric tensors fifj and \ijki- We impose the requirements 

= [/An, 

= ^im ^rnjkl + Tj m \imkl + T£ m \ij m l + T\ m ^ijkm , (14) 

in order to make V(<$>) gauge invariant under G. The first equation implies that the 
matrix fi 2 is proportional to the unit matrix on each subspace irreducible under G. 
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To break the GUT symmetry, there has to be one G-irreducible Higgs representation con- 
tained in <3>j that develops a vev. In order to treat the symmetry breaking appropriately, we 
define the projector H n on this particular G-irreducible subspace (clearly, [T Q ,I1 W ] = 0). 
This subspace is further divided into the subspace of Goldstone bosons and the subspace of 
physical Higgs bosons with projectors P g and P n , respectively ([T a , P g ] = = [f a , P H ]). 
The physical Higgs bosons receive masses of order Mq from V($), the Goldstone bosons 
do not. Using the Goldstone boson mass matrix from eq. (j!2p . the projector on the space 
of Goldstone bosons can be explicitly written down as [ID] : 

(P% = (H* - P% 3 = g*fiv k (J^ vflfj . (15) 

Now we can parametrize the scalars in the following way: 

= v i + & i = v i + H i + G i + Si (16) 

where v, H and G live only on the subspace defined by U n . S parametrizes all the other 
scalar^: 

Si , 

Hi + Gi , 
Hi , 

Gi. (17) 

First let us focus on the subspace that U n projects on. In order to develop a vev on this 
subspace, the parameter defined by U n fi 2 = /i^fl, has to be positive. If this is the 
case, it is convenient to parametrize this part of the scalar potential in terms of physical 
parameters as the Higgs mass, the heavy gauge boson mass, the gauge coupling g and the 
tadpole instead of the unphysical couplings ^ and 

Kjkl = ^i'j'kn'^i^j^'k^-l'l ■ ( 18 ) 

In principle, this is analogous to what is usually done for the SM Higgs potential |JT]. 
Here, however, it is more involved due to the appearance of the general invariant tensor 
Xlj kl . Using essentially eq. ( 1T41) and Schur's lemma, it is possible to rewrite the up to 
quadratic terms of the potential in terms of new parameters Mjj (diagonal Higgs mass 
matrix) and t (tadpole). For the trilinear and quartic terms this does not seem to be 
possible at the level of the Lagrangian. Thus, for the moment, we leave those terms 
expressed by the old parameters A^ fc/ . They have to be eliminated in favor of M^,g and 
M| at diagram level to make our choice of parameters consistent. 

3 To see that the number of Goldstone bosons is equal to the number of broken generators it is also 
possible to define the Goldstone field as G A := g (ife) Vi ^i^^j 



(1-It«)^ = 

= 

U 3 
*J 3 

p9.g>'. = 

13 3 
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Now, including also the scalars on the subspace defined by 1 — U n , which is straightfor- 
ward, the scalar potential can be parametrized as follows 



V($) = tViHi + ^M^^jHiHi + ^tHiHi + ^tGiGi + ^MDijSiSj 



+ 2 V i>H jk iHjG k Gi + -Vi\y kl HjH k Gi + -v i Xf jkl H j H k H l 

+ -ViXijkiHjHkSi + ViXijkiHjGtSi + -ViXij k iGjG k Si 

1 1 1 

+ -ViXijkiHjSkSi + -ViXijkiGjSkSi + -ViXijkiSjSkSi 



+ -^fjki G iGjG k Gi + -Xf jkl GiGjG k Hi + -Xf jkl GiHjH k Hi 



+ -^jki G i G i H kHi + —X™ kl HiHjH k Hi 

1 1 1 

+ -Xij k iHiHjH k Si + -Xij k iHiHjG k Si + -Xij k iHiGjG k Si 

+ -^Xi,j k iGiGjG k Si + -Xij k iHiHjS k Si + -Xij k [HiGjS k Si 

+ -Xij k iGiGjS k Si + -Xij k iHiSjS k Si + -Xij k iGiSjS k Si 
1 

+ —Xij k iSiSjS k Si (19) 



where 



t = -/4 + Xf jkl ViVjV k vi , (20) 

(M^)ij- = -AgfcjUfcU, - ^ X Mmn.VkVlV m V n Ilg , (21) 

= ^gfc» - (A* a (* - n *)k ■ ( 22 ) 

Due to gauge invariance under the SM group (eqs. (J4]) and (|T4j) ), both mass matrices 
are diagonal and proportional to the unit matrix on each SM-irreducible subspace. M\ 
has only non-zero entries on the subspace defined by P n and Mj only on the subspace 
defined by (11 - U H ) (we have defined Xf jkl v k v t = (1 - 1^)^(1 - n w ) ij /A^, fe ^ fc ^ ). 

It is important to see that Mj must have only positive or zero entries. If there are negative 
entries, some of the Si would develop a vev and our formalism would not apply. Strictly 
speaking, we have (Mj)jj < for the SM Higgs doublet that is contained in Si, which 
would exclude it from our treatment. But since in that case the scales involved have the 
strong hierarchy 0(M W ) <C 0(M GUT ), we can safely set the entry to zero here. To do 
this, some of the Xfj kl v k V[ must be fine-tuned against the corresponding (/i 2 (ll — IT^))^- 
in eq. (122]) which is known as the doublet triplet splitting problem, inherent to generic 
GUTs. 



4 For more details of how this reparametrization is done, please refer to appendix! 
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Note that the classical minimum of the GUT-breaking Higgs potential is defined by the 
equation t = 0. However, if we compute higher order corrections, the parameter t = — 8t, 
where St = 0(a) is a counterterm, has to be adjusted in such a way that the renormalized 
Higgs one point function is zero at all orders of perturbation theory. 

The last term in eq. (CQ) to be specified is Cy- The most general Yukawa interaction of 
the chiral Dirac fermion multiplet \l/ with the real scalar multiplet $ can be written as 
follows: 

C Y = -~ (Y^JCV^k + Y^fC^ k ) . (23) 

Here Y k is a complex, symmetric matrix and C = 27270 denotes the Dirac charge con- 
jugation matrix. \& c = is the charge conjugated Dirac spinor (the T refers only to 
Lorentz space). We take \& = 75)^ to be left-handed so that \I/ C will be right-handed. 
Furthermore, due to gauge invariance Y£j satisfies the following relation: 

fl = Y k T a 4- Y k T a 4- Y m T a (2<t) 

mj mi ' im mj ' ij mk ■ \ ) 

In the following we will need to distinguish between fields with the mass of O(Mq) and 
massless fields. We will follow the convention of appendix [A] and use the projectors Pf 
for heavy fields and pf for the light fields. 



2.2 Renormalization 

In order to do a two-loop calculation of the matching corrections, a one-loop renormaliza- 
tion program has to be carried out for the theory. The counterterms are adjusted in such 
a way that all the one-loop Green's functions of the theory are finite. For convenience 
we use the on-shell scheme for the mass parameters of the theory and MS for the gauge 
couplings, the gauge parameters and the fields. The renormalized Lagrangian is obtained 
from eq. ([!]) by the following replacements: 



A a ; ^Z 3i A* , 
pf pf # , 



P?H^^/Z m P?H, 
pfS^^pfS, 
Ml -+ Z 2 Mxi Ml , 
M Fi -> Z Mpi M Fi , 



Pf^ -»■ \/Z& Pf * 




Pi G — > \fZ~Gi P { G , 

PfS^y/ZHPfS, 

Ml -> Zl m M\ , 



6i £ • 



(25) 
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Again, we have used the sub-index i to take care of the fact that there might be several 
SM-irreducible representations for a field that all renormalize differently. No summation 
is performed over that index. Pf and pf are projectors on the various SM-irreducible 
subspaces of heavy and light fields, respectively. Mp is the fermion mass matrix that can 
arise from the Yukawa interactions eq. ( 1231) . Presently, heavy fermions are not included 
in the calculation, so the corresponding renormalization constants are only defined for 
future convenience. 

In the following we list the counterterm Feynman rules that are important for our cal- 
culation. They are obtained by inserting the renormalization prescriptions from eq. (1251) 
into eq. ([T|) and considering the up to quadratic terms. For each counterterm we give an 
expression that is valid to arbitrary loop order in the first line and in the second line a 
more convenient expression that is valid only for one-loop renormalization. We use the 
notation Zj = 1 — 5Zi and t = — St where SZi and St are of order a. All the parameters 
that appear in the equations are renormalized ones. 

Heavy gauge boson: 



Aa Bv 



Heavy ghost: 
^k~ 

Goldstone boson: 

i X 3 

Physical Higgs boson: 



i$AB 

%S A B 
i$AB 



^-^ + \-l)k,K-{Zi-l)eg,MZfZl^-l)Mlg, u 
-5Z^K-(Ml-k 2 )5Z^g^-25Z Mx Ml giMU 



(z£-i)k 2 -(zx^^z* Mx -i)ZMl 

5Z^iMl-k 2 )+{Uz^+Uz^+25Z M ^Ml 



'Z G -l)k 2 -{Z G Z i2 Z^-l)iMl - 1 



5Z G {iMi-k l )+{5Z^+28Z Mx )£Mi + St 



% X 3 



(Z H -l)k 2 -(Z H Z 2 M -l)M 2 H - 1 



8Z H {M 2 H -k 2 )+28Z MH M 2 H + St 
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Mixed counterterm 



i 



W \ t - 1 V k T ki k , 




k 



W 7^% ~ 5Z ^ v ^ki k , 



Heavy fermion: 



(Z%Z Mf -1)M f p l 



I 



k 



j 



iSij 5Z%(M F 



ft)+8Z MF M F P L 



i 



Higgs tadpole: 

X 



iViSt 



In order to avoid clutter with the notation, we have omitted the sub-index i here. From 
the context it is always unambiguous that the SM-irreducible representation of the field 
under consideration is meant. For our calculation we need the renormalization constants 
of all mass and gauge parameters at one-loop. As can be seen from the above Feynman 
rules, the set of equations that is used to determine SZ^, 5Z& and St is overconstrained 
(cf. also refs. (I2J S3]). This provides a useful check for our calculation: we computed 
SZfr — SZ^ 2 from the pole of a combination of the heavy gauge boson and heavy ghost 
propagator as well as from the mixed Goldstone boson heavy gauge boson propagator. 
Both calculations lead to the same result. In the same way St was computed from the pole 
of the physical Higgs tadpole as well as from the Goldstone propagator yielding the same 
result. Furthermore, the MS renormalization constant for the gauge coupling Z g will 
be needed to two-loop order (cf. subsection 12.31) . including corrections from the Yukawa 
couplings. We did the calculation in our framework and found agreement with the known 
result in the literature [2T | l ^ ll5]. 

2.3 Decoupling of heavy particles 

When studying gauge coupling unification, it is most convenient to use RGEs defined 
in mass-independent renormalization schemes, such as MS or DR. In these "unphysical" 
schemes the computation of (3 functions is simplified significantly, which makes them most 
suitable for this application. However, these schemes have the well known drawback that 
the decoupling theorem jl6] does not hold in its naive form. As a consequence observables 
of low energy processes will depend logarithmically on all the heavy particle masses of 
the GUT. This is unacceptable, since it would spoil perturbation theory by the presence 
of large logarithms \jo.(Mq / fi) , where \i is the typical energy scale of the process. The way 
out of this dilemma is to use an effective theory [21H321I1Z], where the heavy particles are 

5 At the two-loop order also more complicated functions of the heavy masses can appear. 
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integrated out at the GUT scale. This means that the dynamical degrees of freedom of the 
heavy particles are removed, which manifestly leads to power-suppressed contributions of 
0(1 /Mq) in the effective Lagrangian. Moreover, the effects of the heavy particles are 
encoded in a multiplicative redefinition of all the masses and couplings of the theory. For 
the case of the gauge coupling this so-called decoupling relation reads: 

ai(yUGUT) = Ca l (A i GUT,«(/iGUT),M h ) a(/i G UT) , 

<* = <* = f^ ' = 1,2,3. (26) 



Here ctj and a stands for the MS gauge couplings in the effective theory (the SM or 
the MSSM) and full theory (GUT), respectively, //qut is the unphysical scale, at which 
the decoupling is performed. At sufficiently high loop order predictions of physical ob- 
servables must not depend on fiQVT anymore. The remaining dependence on this scale 
gives us an estimation of the theory uncertainty of the prediction. ( ai is the so-called 
matching coefficient that depends on all the mass parameters of the particles that have 
been integrated out. They are abbreviated by Mh in eq. ( l26l) . The construction of the 
effective Lagrangian is described in detail e.g. in refs. [181H9] for the case of QCD. It is 
applicable to our setup without modifications. Here we only list the formulas that are 
relevant for the computation of ( ai . For our calculation we find it most convenient to use 
the three-point Green's function with light ghosts and light gauge bosons as external par- 
ticles. Applying Slavnov- Taylor identities to this vertex, the formula for the MS matching 
coefficient reads: 



<- = 1 ¥A) ' i=i ' 2 - 3 - (27) 

The MS renormalization constants for the gauge coupling in the full and effective theory 
are denoted by Z g and Z 9i , respectively. The bare matching coefficients for the light 
ghost-gauge-boson vertex, the light ghost field and the light gauge field respectively are 
given by: 

a = i+ r i ft t c ,(o>o), 

S = i + n°f(o), 

C 3 ° = i + njfto) (28) 

i.e. they are computed from the "hard part"0 of the one particle irreducible Green's 
functions with zero external momentum. All the masses of SM particles are set to zero 
and therefore only masses of O(Mq) appear in the diagrams. The external particles for 
these Green's functions are: two light ghosts with one light gauge boson, two light ghosts 
and two light gauge bosons, respectively. In all cases the relevant group structure has to 



6 For simplicity we will only speak about MS parameters from now on. If a SUSY GUT is considered, 
all MS parameters will be replaced by DR parameters. 



7 This denotes all the diagrams that contain at least one heavy particle. 
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be projected out and in the case of the light gauge boson propagator only the transverse 
part is needed. The index % takes care of the fact that the SM gauge group is not simple 
and labels whether an external gauge bosons and ghosts belonging to 17(1), 577(2) or 
SU (3) have to be taken. 

At this point a remark concerning the /xqut dependence in different renormalization 
schemes might be helpful. In mass- dependent renormalization schemes the threshold 
effect around the GUT scale is obtained correctly by including the effect of the heavy 
GUT particles in the f3 function. As there is no explicit matching performed in such 
schemes, there is no dependence on the unphysical scale /igut at all. Therefore, the ques- 
tion might arise how the /xgut dependence in mass-independent renormalization schemes, 
as in eq. ( 12"T|) . comes about. The answer is that by naively changing the renormalization 
scheme from mass- dependent to mass- independent, low energy observables obtain a log- 
arithmic dependence on an unphysical renormalization scale fi that has not been present 
before. These are the large logarithms log that would spoil perturbation theory in 
the naive mass-independent renormalization scheme. In our (effective theory) approach, 
however, we absorb these logarithms in a redefinition of the gauge coupling and choose 
the unphysical scale /i = /xgut in the vicinity of the heavy GUT masses, which introduces 
an explicit dependence of ( ai on /xgut- 

Before we come to the technical details of the calculation, it is in order to describe all the 
assumptions that have been made about the underlying GUT model: 

• There is no trilinear scalar coupling in V(Q). 

• There is only one vev of O(Mq) in the theory. 

• There are no heavy fermions in the theory. 

• The heavy gauge bosons decompose in SM-irreducible representations with a com- 
mon mass. 

• The GUT-breaking Higgs decomposes into three SM-irreducible representations (+ 
Goldstone bosons) at most. They can all have different masses. 

• The other scalars in the theory decompose in SM-irreducible representations that 
have a common mass (+ light scalars). 

• The light particles in the theory can decompose in arbitrarily many SM-irreducible 
representations. 

As can be seen, the main limitation comes from the number of heavy degrees of freedom 
in the theory. The above constraints are designed such that the resulting formula for £ a . 
is applicable to the simplest GUT, the Georgi-Glashow model, yet keeping the calculation 
as simple as possible. The computational framework for our calculation is set up in such 
a way that it can be generalized to more heavy degrees of freedom, in order to apply it 
to SUSY GUTs in the future. 
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Given the large number of Feynman diagrams, an automated computation is indispens- 
able. The diagrams were generated with QGRAF [50] and further processed with q2e and 
exp [5TTI52] . In the next step we used a FORM [53] implementation of the two-loop topolo- 
gies of ref. [51] by the authors of ref. [13] and also the FORM packages MINCER [55] and 
MAT AD [56]. Let us emphasize that no assumptions about the mass hierarchies of the 
heavy particles have been made. Therefore, the result is valid for arbitrary numerical 
values of the mass parameters as long as their mass splitting is not to large which would 
lead to power enhanced contributions and spoil perturbation theory. The reduction of the 
group theory factors, which was the most time consuming part, has been automated and 
implemented in FORM. The relevant group theory and notational conventions are collected 
in appendix [A] of this paper. 

Although the present calculation is just general enough to be applied to the simplest 
GUT, the number of Feynman diagrams for the two-loop Green's functions described 
above already is considerable. For the light gauge boson two-point function it amounts 
to 6278, whereas for the ghost-gauge-boson vertex and the light ghost two-point function 
we have 4109 and 374 diagrams, respectively. Sample diagramazl for all three processes 
are depicted in fig. CD 

In order to subtract the sub divergencies that occur in these two-loop Green's functions, the 
mass and gauge parameters as well as the gauge coupling that appear in the corresponding 
one-loop Green's functions have to be renormalized as described in subsection 12.21 We 
have performed our calculation for arbitrary gauge parameters and verified that they 
cancel out in the final result, which is a powerful and highly non-trivial check of the 
calculation. The result for ( ai is available in general form, i.e. with group theory factors 
and couplings not specified to a particular Lie group or model. In the next section we 
will assign definite values to these quantities in order to show the application of the result 
exemplarily. Note also that we always need three different sets of group theory factors for 
i = 1, 2, 3, respectively. 



3 Numerical study for the Georgi-Glashow model 

In the following we perform an academic study of our results applied to the simplest 
possible GUT, the Georgi-Glashow model [54] . Although this theory is ruled out experi- 
mentally [2] , we use it as a toy model to demonstrate some simple numerics. It is based 
on the gauge group SU(5) and the SM fermions sit in the representations 

5 (3,l,|)©(l,2,-i), (29) 
10=[5x5] a (3,l,-|)©(3,2,|) a ©(l,l,l). (30) 

8 The figure has been created with help of the IATgX package AXODRAW [57] 
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Figure 1: Sample two- loop diagrams that appear in the calculation of ( ai . The first 
line shows the process Af* —> A b j contributing to 11^(0). The second and third line 

depict c a% -)> c h and c a ' ->■ c h + A% contributing to U° c f (0) and r^ tci (0,0), respectively. 
Colored (bold) lines represent fields with mass of O(Mq) and black (thin) lines massless 
fields. Furthermore, curly lines denote gauge bosons, dotted lines ghosts, dashed lines 
scalar fields and solid lines fermions. Goldstone bosons are marked green (light gray, 
short-dashed) , physical Higgs bosons red (gray, long-dashed) and other heavy scalars blue 
(dark, short-dashed). Note also that two identical lines in one diagram need not have the 
same mass because of the non-degenerate mass spectrum. 



where we have displayed their decomposition into SM multipletsjj The gauge bosons as 
well as the S*?7(5)-breaking Higgs bosons E live in the real 24 representation: 

24 (8, 1, 0) © (1, 3, 0) © (1, 1, 0) © (3, 2, -f ) © (3, 2, §) . (31) 

In the case of the gauge bosons the first three parts of the decomposition constitute the 
light SM gauge bosons and the last two parts the heavy gauge bosons with a common mass 
Mx- For the GUT-breaking Higgs bosons the first three parts represent the three physical 
Higgs bosons with masses Ms, 2 Ms and M24, respectively. The last two multiplets give 
the Goldstone bosons with the unphysical mass Mc (cf. eq. (fl2l)). Note that we write 
the field E as a 24-dimensional vector multiplet in order to be consistent with our notation 

9 The first and second number label the SU(3) and SU(2) representations, respectively and the third 
number is the hypercharge of the multiplet. The index a indicates that the antisymmetric part of the 
representation has to be taken. 
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in section [2J and not as a hermitian 5x5 matrix as usually done. Therefore, the vev of 
E simply is (£) = v, where v is a 24-dimensional vector with a single non-zero entry in 
the 24th component. Additionally we have a complex scalar in the 5 representation that 
decomposes according to eq. (129j) and contains the SM Higgs doublet as well as a heavy 
colored triplet with the mass Mh c . 

For the convenience of the reader it might be helpful to explicitly give the parametrization 
of the quartic scalar coupling A^m from eq. (fl3l) for the Georgi-Glashow model. It splits 
up into three parts. The first part is a quartic coupling of the 24 Higgs, the second one 
a quartic coupling of the 5 Higgs and the last one is a mixed 5 — 24 coupling [40l[58] : 

X 2 a % 5 = A sTr(T a T^T 5 ) + l -B {8 aP 8 lS + 5 ai 5 0S + 5 aS 5 M ) , 

>Hjki = 3 b + 8 ik Sji + S a S jk ) , k,l = 1, 10, 
^J! = cCrV + rVOy, a,/3, 7 ,5 = l,...,24. (32) 

We have used the symmetrized trace 

sTr(T Ql ...T°") = ^^Tr(T Q -w...T°-(")) (33) 

7T 

where the sum is over all the permutations of the indices. Furthermore, it is noticeable 
that the indices i,j,k,l run from 1 to 10 although they belong to the fundamental 5 
representation. This is because the corresponding scalar is complex and we have written 
it as twice as many real scalars that are transformed by the 10 x 10 generator matrices r: 



i Im(T a ) i Re(T Q ) 
-i Re(T°) i lm(T a ) 



(34) 



where in this section T a is the 5x5 generator matrix in the fundamental representation. 
Inserting eq. (|32|) into eqs. (l2Tj) and ( )22|) . we obtain the scalar mass matrices. Addition- 
ally we need to impose the tree-level fine-tuning condition /i| = t^cv 2 , where /i| is the 
quardatic term of the 5 Higgs in eq. ( ITB"]) . in order to obtain massless Higgs doublets. Note 
that in principle one would need to calculate the one-loop fine-tuning condition in order 
to obtain massless Higgs doublets in a two-loop calculation. However, the light Higgs 
doublets show up at the first time in the two-loop Green's functions in the matching 
calculation so that it is sufficient to use the tree-level fine-tuning condition. This gives us 
the relations between the physical scalar masses and the parameters in eq. ( 132]) : 

Ml = ^Av\ M 2 24 = l( J L.A + B)v\ = i c v 2 . (35) 

The vev is connected to the physical gauge boson mass by 

Ml = ^g 2 v 2 . (36) 
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The Yukawa interactions of the Georgi-Glashow model [3UI10] are obtained by inserting 
the Yukawa matrix 

yn _ f —YlJ^ijklm T^T^S^n 2z Yjj Im^^jSl 

- V 2iY£lm(TS)S ln 

into the general Yukawa Lagrangian eq. fl23l) . Here s = (I,s) and r = (J,r) are multi- 
indices, where /, J stand for the generation indices of the SU (5) Yukawa matrices Y u 
and Y D . The indices s, r = 1,...,29 run over {a,j} and {/3,k}, respectively. Note that 
we have written the fermions of the 10 representation as a 24-dimensional vector instead 
of an antisymmetric 5x5 matrix as usually. The Clebsch-Gordan coefficients for this 
transformation are given by the following equations: 

10° = V2 Wij , lOij = -y/2 % Im(7£) 10° , (38) 

where 10jj is the usual antisymmetric 5x5 matrix with the normalization as in ref. [34] . 
Furthermore, S = -^(1, zl) is a 5 x 10 matrix and e^kim is the totally antisymmetric 
tensor with €12345 = 1. As can be seen from eq. (1371) . the chiral fermion multiplet \P 
from subsection 12.11 is written as a 3 (24 + 5) = 87-dimensional vector for the case of 
Georgi-Glashow SU(5) model. 

Using the definitions from this section, we computed the numerical values of all the group 
theory factors that appear in our general result. Furthermore, we set Vckm = H and 
kept only the third generation Yukawa couplings yt and y^. We obtained three two-loop 
formulas for ( ai (i = 1, 2, 3) that depend on the parameters 

a(fi G u T ), y t {fJ>GVT), 2/&(a*gut), M x , M Hc , M s , M 24 , a*gut ■ (39) 

The Mathematica package that contains the expressions can be downloaded from 

http : //www-ttp . particle . uni-karlsruhe . de/Progdata/ttplO/ttp 10-46/ 

In order to examine the numerical impact of the two-loop matching corrections in this 
model, we have implemented a RGE analysis in Mathematica. Since in this model the 
gauge couplings do not unify, we just focus on examining the reduction of the decoupling 
scale dependence, as an illustration of our results. We start with the precise values of 
the three gauge couplings at the electroweak scale. They are obtained from the effective 
weak mixing angle in the MS scheme [59] , the QED coupling constant at zero momentum 
transfer and its hadronic [60] contribution in order to obtain its counterpart at the Z-boson 
scale, and the strong coupling constant [61] These quantities need to be transformed to 
a six-flavor theory, which is described in detail in ref. [H]. Our starting values are then: 

1/(128.129 ±0.021) , 
0.23138 ±0.00014, 

0.1173 ±0.0020. (40) 

10 Wc adopt the central value from ref. [61], however, use as our default choice for the uncertainty 0.0020 
instead of 0.0007. 




sin 2 ( 6)i MS (Mz) = 

af\M z ) = 
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These quantities are related to the three gauge couplings via 




3 cos 2 e( 6 )> MS 



sin 2 e( 6 )< MS 





(41) 



which holds for any renormalization scale /i. We also need the W and Z boson pole 
masses My/ and Mz, the top quark and tau lepton pole masses M t and M T and the 
running bottom quark mass mj^ s . For the convenience of the reader we also specify their 
numerical values [59 | l62 | l63]: 



The corresponding uncertainties are not important for our analysis. These parameters are 
converted to six-flavor theory using RunDec [M] and then used to compute the starting 
values for y t , yb and y T at the electroweak scale. The RGE running in the SM was imple- 
mented at two loops [4^114311651166] for the electroweak sector and at three loops [671168] for 
QCD. We take into account the tau, bottom and top Yukawa couplings and thus solve the 
coupled system of six differential equations. Since the quartic SM Higgs coupling b enters 
the equations of the Yukawa couplings starting from two-loop order only, we neglect its 
contribution. After taking into account the two-loop decoupling relations, we compute 
the running from /iqut to the Planck scale using three-loop RGEs for the gauge coupling 
and one-loop RGEs for the Yukawa couplings. The RGEs are obtained by inserting the 
general expressions for the Yukawa and scalar couplings (eqs. (137|) and (1321 ) as well as 
the numerical values for the group theory factors into the general formulas of refs. [21~||4"5] 
(see appendix ICl for the details). 

In figure |2] the dependence on the decoupling scale of a(10 18 GeV) is shown. Since only 
for QCD the full three-loop function could be implemented and there is no unification 
of gauge couplings anyway, we took o;(/^gut) = CmH^gut) «3(/^gut) as a starting value 
for the gauge coupling above the GUT scale. For illustration we use the following set of 
mass parameters: 




M w 
M z 
M t 



80.398 GeV, 
91.1876 GeV, 
173.3 GeV, 
1.77684 GeV, 
4.163 GeV. 



(42) 



M x 

M E 
M 24 



10 15 GeV, 
4- 10 13 GeV, 



10 14 GeV , 
6 • 10 13 GeV 



(43) 
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0.0178 
0.0177 
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¥ 0.0174 
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0.0172 
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log 10 (y"GUT/GeV) 

Figure 2: Dependence of a(10 18 GeV) on the decoupling scale //gut- The red (dotted), 
green (dashed) and blue (solid) lines correspond to the one-, two- and three-loop analysis, 
respectively. For the three-loop curve also the experimental error band with 5a s = 0.0020 
has been indicated. 



which are chosen to obey the restriction Mx > for % = H c , S, 24 . Otherwise the scalar 
self-couplings easily become non-perturbative and blow up the gauge coupling above the 
GUT scale. The scale dependence is shown for n-loop running and (n — l)-loop decoupling 
with n = 1,2, 3. We observe a dramatic improvement when going from n = 2 to n = 3. 
In particular the three-loop corrections can be larger than the experimental error band 
depending on /xqut- Note also that for n = 2 choosing /^gut naively as a mean value of 
the GUT masses which would be of C(10 14 GeV) in our case is not a good choice. The 
described qualitative behavior does not depend much on our choice of the GUT masses. 
Though the numerical effect of the two-loop matching is already significant in the Georgi- 
Glashow model, we emphasize that in certain models that contain large representations, as 
e.g. the Missing Doublet Model [TSIITB]. we expect these corrections to be even larger [T4"j . 
Our goal for the future, of course, is to generalize the formula for ( ai to make it applicable 
to these models. 

To provide a check of the result for the Georgi-Glashow model, we have verified ana- 
lytically that the matching coefficients Ca^cirr) exhibit the correct //gut dependence. 
This can be derived from the knowledge of the two-loop j3 functions of the SM and the 
577(5) model by computing the derivative w.r.t. £gut = hi^GUT) of eq. (12"E|) . Solving 
the resulting differential equation order by order, we arrive at a general formula0 for the 

n For simplicity we neglect the Yukawa corrections in this formula. However, the generalization is 
straightforward. Of course, in our analytical check we took care of them too. 
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/U-gut -dependent terms in Cai(/-*Gur) 



C^gut) = 1 + ^ GUTj - A,) *gut - C (M h ) (44) 

7T L J 

2 

- /3o) 2 4ut + UPi- - ^(Mh) - /3 ) t GUT + Ci(M h ) 



+ 



«Qgut) Y 

7T J 



Co and Ci are /XQUT-independent terms that depend only on the heavy GUT masses. The 
j3 function coefficients are defined by: 

fc=0 

and similarly for We find agreement in the /igut dependence of our explicit calculation 
with the form of eq. 



4 Conclusions and Outlook 

As experimental accuracy for a s , a em and sinB is increasing, also theoretical unification 
analyses must improve their precision in order to find better exclusion limits for GUTs. 
Therefore, we have performed a first step towards the calculation of the two-loop matching 
corrections for the gauge couplings at the GUT scale in a framework that aims at making 
as few assumptions on the underlying GUT as possible. The assumptions that were made 
can be found in subsection 12.31 The result is general enough to be applied to the simplest 
GUT, the Georgi-Glashow SU(5). The numerical impact in this model was found to be 
larger than the current experimental uncertainty on a s depending on the choice of the 
decoupling scale //gut- We expect even larger effects in GUT models that contain large 
representations, as the so-called Missing Doublet Model. Moreover, the two-loop matching 
coefficients provide a significant stabilization of our predictions w.r.t the variation of the 
decoupling scale. 

Furthermore, we have described in detail the proper treatment of the gauge fixing and 
the renormalization procedure for this calculation. In this context also the issue of Higgs 
tadpoles in theories with spontaneous symmetry breaking was discussed in a general man- 
ner. The appendix contains an in-depth introduction of the group theoretical framework 
that we used. There we also give many useful reduction identities that are essential for 
performing multi-loop calculations in GUTs and can also be applied to SUSY GUTs. 

We consider the result computed in this paper only as an intermediate step towards a 
more general calculation of ( av Particularly, we are interested in performing a consistent 
three-loop RGE analysis for SUSY GUTs, where two-loop matching at the GUT scale 
is needed. In order to apply our result to the simplest SUSY GUT, the minimal SUSY 
SU(5), several generalizations are needed: 
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• Add a trilinear term to the scalar potential eq. (|T3|) . 

• Add two massive Dirac and three massive Majorana fermion SM-irreducible repre- 
sentations as well as all possible kinds of Yukawa interactions for them. 

• Increase the number of SM-irreducible representations within TZ n by one in order 
to cover the case of the non-renormalizable version of minimal SUSY SU (5) 

• Allow for unitary mixing matrices in the interactions of heavy Dirac fermions with 
gauge bosons. 

• Convert the result to the DR renormalization scheme. 

Though there is quite some work to be done, to achieve this, we have already built a solid 
basis to start with and have overcome many of the main difficulties that were expected 
to occur. 
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Appendix 

A Group theoretical framework 

When performing multi-loop calculations for GUTs, one necessarily has to deal with the 
reduction of the color structure^ [ZDl23j|IIJ|l5j|69j[70] ' ^ e a i m * s t° reduce all the color 
factors of the diagrams to a set of basic invariants. In this way nontrivial cancellations 
among different diagrams are possible without inserting the actual numerical values for 

12 For convenience we will use the terms "color structure", "color factor" etc. in this appendix following 
QCD terminology. However, obviously our formalism is not restricted to SU (3), but is meant to be applied 
to GUT groups. 
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those invariants. In this appendix we develop a notational framework that is appropriate 
for this task and also give some useful reduction identities that will hopefully prove helpful 
also for future calculations. For similar reduction algorithms and identities applied to 
unbroken gauge theories see e.g. refs. [HIES] and references therein. Here we focus on 
spontaneously broken gauge theories. 

We start with a simple GUT group G that is broken to a (in general not simple) gauge 
group Ylk Gfc. In the following the G& are called group factors. The generators of G in a 
general reducible representation 1Z are denoted by T a in this appendix^!. They fulfill the 
commutation relation 

[T a 1 T p ]=if a ^T rf , (46) 

where f al3/y are the structure constants of G. In order to distinguish between broken and 
unbroken generators, we use the notation of subsection 12.11 

{a} = J2{Ai} + = {A} + {a} (47) 

i i 

where A\ label the broken generators of G belonging to the Yl k G^-irreducible subspace 
labeled by i. If there is only one TJ fc Gfc-irreducible subspace in the adjoint representation 
of G, we can omit the sub-index % in Aj. In contrast, a, label the unbroken generators 
belonging to the subgroup G,. As the Gj are regular subgroups of G, also the following 
cummutation relations hold: 

Jjia^ jAjj _ ^ jaibjC k rpc k ^g-j 

where f aih i Ck = unless i = j = k. Furthermore, because the subgroup Gj is closed, 
we have f aib i Ak = for all i,j,k. Otherwise the commutator [T a %T fej ] would contain 
terms proportional to the broken generators T Ak . Note that if not explicitly stated, 
the sub-indices of the indices a,i,bj..., Ai, Bj... are not summed. A repeated index 
a,b... or A,B... without sub- index means that the sub- index has been summed over. The 
representation TZ that T a is defined on generally is reducible under G. It decomposes into 
G-irreducible representations where the gauge bosons, fermions and scalars of the theory 
live in: 

n -> © nX = nA ® nH ® nSl ® nSu ® - ® ® nTl1 ® - • ( 49 ) 

IZ^ stands for the adjoint representation and 7Z n for the representation of the GUT- 
breaking scalar. The other symbols represent the irreducible representations for the scalars 
and fermions, respectively, numbered by roman numerals for convenience. We define 
projectors on these subspaces denoted by IP with x = A,H,Si,Sn...,J-'i,J-'ii.... Clearly, 
[IP,T a ] = holds for all x and a. Each of those representations decomposes further 
under f] fc G fc : 

nX ® n n = ni®K x 2 @Kl ©... (50) 



13 In the main text we employed the symbol T a only for the fermion representation. In this appendix, 
however, we will use the symbol for a generic generator of the gauge group. 
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with projectors P x and p x . We use a capital P to denote that the respective Ylk G*r 
irreducible representation contains fields with a mass of O(Mq) and a lowercase p for 
projectors on a subspace with massless fields. Because the Lagrangian is still invariant 
under Ylk Gfc after symmetry breaking, each of those Ylk G^-irreducible subspaces can 
be assigned to a definite mass of the respective field. The indices n = 1,2,3... label the 
Gfc-irreducible representation in TZ X . The projectors fulfill [P x ,T ai ] = = [p^,T ai ] for 
all x,a,n and i. Furthermore, we have 

E^ = I>n + E P n =P X + P X = TI X (51) 

n n n 

where g can be P or p depending on n. Armed with these definitions, we can define 
some basic Dynkin indices I 2 {---) and Casimir invariants C 2 {...) which have real numerical 
values for a given group. 



ra/3 



Tr(U x T a T p ) = I 2 (U X ) 8 C 



■ 



Tr(g x n T ai T bi ) = I 2 {g x n ) 1 5 a 

Tr(g x T A >g x m T B >) = I a (£, g x J S A ^ 

U x T a T a = C 2 (T1 X ) IF , 

g x T ai T ai = C 2 (g x Jg x n , 

g x n T A >g x m T A * = C 2 (g x n ,g x m Y g x n . (52) 

Again g can stand for either P or p. Furthermore, we can define the dimensions of the 
various irreducible representations by 

A x = Tr(ir), 
K = Tr(P x ), 

d x = Tr(p^), (53) 
and specifically for the adjoint representation: 



A- 4 EE 5 



a a 



^)A — £A n A n 

J n — u i 

d A = 5 a n a n _ ^ 



which gives us the relations 



/ 2 (IP) A A = C 2 (IF) A 

i 2 (PZYdf = c 2 {p x yp. 



X 

J 



II ■ 



i 2 ( P x n ydf = c 2 ( P x n yd x n , 

I 2 {P^g x jD A = C 2 {P x n ,g x J D x , 

h{p x n ,Q x jDf = C 2 (p x nJ g x m yd x n . (55) 
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Let us emphasize again that no summation over the sub-indices % that label the \\ k G^- 
irreducible representation is implied here. However, we introduce the convention that 
omitting this sub-index implies summation: 

c 2 (...) = ^c 2 (...r, 

i 

m 

i 2 (... g *...y = y.^a-lC-y- (se) 

m 

Note also that the invariants I 2 (H X ) and C 2 (U X ) have only been defined for convenience. 
Actually they can be decomposed into other "more elementary" invariants: 

i 2 (ip) = i 2 (p x y + i 2 (p x ) i = i 2 (p x 1 p x ) i + i 2 (p x 1 p x y + 2i 2 (p x ,p x )\ (57) 

C 2 (IP) = C 2 {p x ) + C 2 (p x ,p x ) + C 2 (p x , P x ) = C 2 {P X ) + C 2 (P x ,p x ) + C 2 (P X , P x ) , 

where the summation convention introduced before has been used. The right hand side 
of these equations does not depend on i anymore due to the summation of the projectors 
over the full representation space TZ X . Since in the calculation of matching coefficients 
one has to distinguish between heavy and light particles, it is more convenient to use the 
"more elementary" invariants on the right hand side. 

The definitions that have been introduced are general enough to be applied to all color 
factors that appear in our calculations. However, some of the fields live in representations 
that deserve special attention. First let us focus on the adjoint representation. The 
generators here are defined by 

(T%) M = {n A T a ) M = -if^ . (58) 

Clearly, the operators P A and pf project on the subspaces with indices Aj and a,, re- 
spectively. Because now x = A in eq. (I55p and / a « 6 J A * = 0, things simplify for the adjoint 
representation. In fact it is sufficient to define four invariants: 

/ 2 (n^) s a/) , 

hipfy ^ , 

h{pfy s aibi , 

h{Pf,pfy 8 A ^ (59) 

where we followed the notation of eq. (152]) . We keep the redundant sub-indices i and j 
in lines 2 and 4 of this equation in order to be consistent with the notation of eq. fl52l) . 
Note that f a i B J c k = q f or j ^ k because the indices j and k are used here to label 
the Ylk Gfc-irreducible representation of the generator TjJ* . There is another quadratic 



ya7o JP'yo 
jdiCidi jbiCidi 

["'■"■ f'' ( P 

J ' '' '' ./'''' ' ' 
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Casimir invariant for the adjoint representation that can be expressed through these and 
one that vanishes: 



A,CD rBjCD 



2i 2 (p A , P A y 



S A i B i 



jAiCD jbjCD _ g 

Furthermore, there are relations between these invariants: 

h{ii A ) = i 2 (p A y + i 2 ( P fy 

At the two-loop level these relations are not sufficient to reduce all adjoint color factors 
because products of up to six structure constants with various contractions can appear in 
the diagrams. Using the Jacobi identity for / Q/?7 , we can derive relations for products of 
three contracted structure constants that have three open indices: 



(60) 



(61) 



jceSe j?/3e(f> j?-y4>8 

jciidiei jbieiSi j(Hfi<k 

jaiDE jhEF jCiFD 
jdiDE jbiEF jCFD 
jaiDE jBjEF jCjFD 

jaidiCi jBjeiF jCjFdi 

j- ai DE jB k Ef 3 jCkfjD 

jAjDe jBjeF jCjFD 

jAjDE jBjEF jCjFD 



-h(H ) f 



A\ fafi^i 



A\i faciei 



Uptyf 



•A\i faibiCi 



h(P A yf 



2 

0. 
1 

2 

\h{pfyr B ^ 



i 2 (u A )-2i 2 (p A P A y p B *°i 



2i 2 (P A , P A ) k -6 l] i 2 (p A y 
i-Aif.i> A yf All ' r . 
h(u A )-3h(P A , P A y}r 



f 



o-iB k C k 



(62) 



These relations are sufficient to do all the reduction for two-point and three-point Green's 
functions for the adjoint representation at the two-loop level. 

Next, let's turn our attention to the representation lZ n where the GUT-breaking scalar 
field H and the Goldstone field G live in. Some peculiarities occur here due to the 
appearance of the vev v. 7Z n decomposes under FJ fc G& into the part of physical Higgs 
fields and a part of Goldstone bosons: 



H 



n\®n G 2 



(63) 



As already explained in subsection 12.11 the explicit form of the projector on the subspace 
is given by [3D]: 



if = gf*v (^) «f* 
\ 1V1 X/ AiB, 



p G = ypf 



(64) 
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where the (diagonal) gauge boson mass matrix 

{Ml) AA = g 2 vf A >f B >v ee MIS** 



(65) 



has been used. The antisymmetric generators in the real representation TZ n have been 
denoted by T a . Accordingly, the projectors on the subspace of physical Higgs bosons can 
be written as 



5>. 
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p n = n rt -p 
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M\j AB 



vT 
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(66) 



where Pf 1 projects on TZf-. With these definitions and using eq. (1481) as well as the 
antisymmetry of T a , we already can derive a useful reduction identity for an invariant 
tensor that appears frequently: 



v f*f B if Ck v = A (Ml - Ml. + Ml )f A ^ c * . 
2g 1 3 k 

One important property follows from T ai v = 0: 

f a *f A >v = [f a %f Aj ]v = -(T%) AjBj f B iv. 



(67) 



(6* 



i.e. T ai acts like the adjoint generator on the subspace of Goldstone bosons. This leads 
to various relations between invariants in 1Z n and in TZ A : 
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(69) 



There are also two nontrivial important reduction identities that involve both types of 
invariants: 
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where x G {H, Si, Sn, Jt_, Jri...}. Additionally, w.l.o.g. we now define Pf 1 to be the 
operator that projects on the subspace where v 7^ (i.e. {P\')ij = -^r is a matrix with a 
single non-zero entry in the component (k, k) where vj~ 7^ 0). Then because of T ai v = 
and vT Ai P^ = vT Ai P 6 P^ = 0, any invariant that contains Pf 1 alone or together with 
some other Pj 1 vanishes: 



H tdH\ 



0. 



(72) 



Since the Higgs mass matrix Mfj (cf. eq. (12 ip ) commutes with all T a \ it is diagonal and 
proportional to the unit matrix on each fj fc G^-irreducible subspace. Therefore, it can be 
written as: 



M 2 H 



(73) 



where M]j. are masses of the physical Higgs bosons and particularly P s M 
the Goldstone theorem. 



due to 



Next we will give some useful reduction identities that involve the quartic scalar coupling 
\lj kl = Xi'fk'i' PftiPfj ^k^vii w hi cn is a totally symmetric invariant tensor under G. 
These identities can be derived by using eq. Q14p and Schur's Lemma as well as the 
definition of the Higgs mass matrix in eq. (12 ip . Some of these identities can be obtained 
by multiplying eq. (fl4|) by and performing derivatives w.r.t $ m . They are used 

to eliminate the coupling X^ kl from the result by expressing it through the physical Higgs 



masses. 
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(74) 



In the same way we obtain some important relations for the coupling A that is not re- 
stricted to the space lZ n : 

X ijkn >(vf A ) iVj (g s ) k , k (g s ) in = (M 2 s f A ) kl - (f A M 2 s ) kl , 
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with x,y, z G {H, Si, Sn, ...}. In the same way we can also derive reduction identities for 
the Yukawa coupling Y n using eq. (T2i 



Ti{g x Y n g y k Y m ^){g z T^T a % m = Tr(g x Y n *g y Y m T a *T a >)(g z ) nm (76) 
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Here x,y G {JF\, JFn, ...} and z G {H, Si, Su, ...}. The list of reduction identities may 
not be exhaustive, but it contains the most important relations. In a similar fashion also 
other identities can be derived using the invariance relations eqs. (1T41) and ( 1241) for the 
invariant tensors. 

In the following we briefly sketch the algorithm that is used for the reduction of all the 
color factors in the diagrams: we have written a FORM program that treats the color factors 
for each individual diagram and reduces them to a basic set of invariants. In a first step 
all reduction identities that involve the quartic scalar coupling Ayjy are applied to a given 
expression. After that any expression will contain traces of strings of generators T ai , 
T Ai , projectors g x n and Yukawa matrices Y n . The adjoint indices and Ai are either 



28 



contracted with each other or with some structure constants. First all the contracted 
adjoint indices are removed, i.e. traces of the form 

Tr(...T a \..T a \..), Tr(...T A \..T A \..), (77) 

by applying the definitions of the quadratic Casimir invariants eq. fl52|) . If the generators 
are not next to each other, we commute them until they are and eventually arrive at 
traces that contain no more contracted indices. Next, expressions of the form 

jCidibi rp/ J^a-i J^bi ^ JCxAiBi rjW <J^ 4 ^Bi \ (78) 

are reduced by using 

jCidibirpdirpbi _ 1 JCidibi jVpaj^&ij _ i ^pA ^ i rpci (79) 

and eq. flTOj) . Again generators that are not next to each other are commuted. Expression 
that contain the projector Pf are treated separately. We insert the explicit form of Pf 
(eq. ( l64|) ) and write the traces in the form v...v, where stands for a string of generators 
T ai , T Ai , and projectors g^. Here we additionally make use of the relations T ai v = and 
Pj^v = (for i ^ 1) in order to eliminate all generators inside the string that have a 
lowercase adjoint index. After that all color factors involving the Yukawa matrix Y n 
are reduced to a basic set of invariants using eq. fl76l) . Finally, we are left with various 
contractions of structure constants which are expressed by the respective invariants using 
eqs. ( 1591) and ( l62i) . The actual program is slightly more complicated than described above 
and one needs to introduce repetitive control structures because not all the reduction can 
be done by a single run. However, the basic procedure is as described. 



B Reparametrization of the scalar potential 

In this appendix give some details on how the parametrization of the up to quadratic 
terms in the scalar potential (eq. (Tl9l) ) arises. We start with eq. (fT3j) and insert the 
decomposition of the scalar field 

$i = Vi+Jh+ + (80) 



where v is the vev with a single non-zero component, H the physical Higgs field, G the 
Goldstone field and S a field, representing all the other scalars, present in the theory. 
For clarity we have also indicated the representations where the individual fields live in. 
Considering only the up to quadratic terms, we arrive at the following expression: 

V{*) = mu + ^MlHiBj + \MlG i G j + ^Mp^ + (9($ 3 ) (81) 
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where we have defined the quantities 

U = 
M 2 = 

Here we already have used the fact that, due to gauge invariance (eq. ( Fl4|) ). the matrix 
^ 2 is proportional to the unit matrix on each G-irreducible subspace and therefore 

l4j v i G j = > PijViSj = , tfjHiGj = , 

tfjHiSj = 0, tfjGiS^O. (83) 

Also due to gauge invariance (eq. ([HJ) ) we have XijkiVjVkVi ~ Vi, because the matrix 
K defined by = \ju v k v i is diagonal and proportional to the unit matrix on each 
SM-irreducible subspace. Hence 

KjkiGiVjV k Vi = = XijkiSiVjVkVi (84) 

since (P G )ijVj = = (P s )ijVj. The parametrization in eq. (18 ip has the disadvantage that 
the masslesness of Goldstone bosons is not manifest there due to the appearance of an 
explicit mass matrix for the Goldstone field. Note furthermore that on the subspace lZ n 
the parameter /i 2 is redundant, because it depends on the tadpole term t, which has been 
chosen as a physical parameter in the Lagrangian. Therefore we want /i 2 also to disappear 
from the mass matrix of the physical Higgs bosons H by trading it for t. To solve these 
issues, we first rewrite the tadpole term by observing that t ~ Vi due to arguments already 
given above. Therefore we define 

ti = tVi, t = — (-rfjViVj + -XijklViVjVkVi) 

= + fo^ X tjkl V i V 3 V kVl ■ (85) 

The classical minimum of the potential can be found by setting t = 0. In a quantum 
theory, however, we must keep t as a counterterm, as explained in subsections 12.11 and 
I2~2l Next we turn to the term \M^ j G i G j in eq. ([SIT) . From the Goldstone theorem 
it follows that this mass term must vanish at the classical level. But since we are also 
interested in quantum corrections, care is needed here. We calculate the term by applying 
the explicit form of the Goldstone projector P s which is given in eq. ( I15D to M 2 . Using 
also the gauge invariance relations in eq. ([14"]) . we can show the identities 

(P G )ij\jklmVlV m = 7 ^{P G )ikXjlmnVjV l V m V n , 

{P%A = ^(P e ), fc /i>^/ = (P e ) lfe /i?,. (86) 
Therefore the following identity holds: 

{P% 3 M% = t{P% k . (87) 



1 



-HijVj + -XijkiVjVkVi , 

+ ^X ijk iv k vi . (82) 
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As expected from the Goldstone theorem, the term vanishes at the classical level, but con- 
tributes as a counterterm to the two-point function of the Goldstone field via t. Similarly, 
we calculate ^MfjHiHj by applying the projector P n = U n — P g to M 2 : 

= t(P% k + M 2 k -t(U n ) ik 

= t(P% k + (M 2 H ) jk (88) 
where we have used eq. ( IHTj) and defined the physical Higgs mass matrix 

(M 2 H h = ^klVkVl-^X^VrnVrnVnU*. (89) 

Note that P s Mjj = and therefore the Higgs mass matrix, as we have defined it, can 
have non-zero entries only on the subspace 1Z . Also the parameter /x^ does not appear 
anymore in the definition of Mfj, as desired. For the fields Si no peculiarities occur, since 
they have no vev. Their mass matrix is essentially given by M 2 : 

(Ml)a = \4ki^vi-^ 2 it-^ H ))ij- (90) 
The up to quadratic terms of the scalar potential can now be written as 

V{$) = tv^ + ^Ml^H^UH^ + hG^^^M 2 ),^ 
which is the appropriate parametrization for renormalizing the theory. 



C Three-loop gauge f3 function for the Georgi- 
Glashow SU(5) model 

For performing a consistent three-loop RGE analysis, apart from the two-loop GUT 
matching corrections also the three-loop gauge (3 function for the Georgi-Glashow model 
is needed. The authors of ref. [2TJ give a general formula for the gauge /3 function of a 
general single gauge coupling theory. Specifiying the group theory factors that appear 
there to the Georgi-Glashow model and insering the scalar self-couplings from eq. (1321) . 
as well as the Yukawa coupling from eq. (1371) into their general result gives us the desired 
(3 function including scalar self-couplings and Yukawa corrections: 



31 



Ida 

YdtAv 



40 / a \ 2 

T V4W 
1184 ( a n3 

v47T, 



(92) 



15 



1007357 / a \4 



9 /^\ 2 
2 Utt/ 



^1 

47T 



47T 



1080 V4tt/ 



+ 



1323 /y f y 3617/y fc \2 
4 UtJ 10 Utt/ 



155 A 11 6 



96 (4vr) 2 20 (4tt) 



51 / 
4 Utt 

493 



A 2 47 AB 



125 £ 
12 (4tt) 2 " 

839 y 2 y 2 
8 (4tt) 4 

1 6 2 



25 c 
T(4tt) 



65 5 2 851 



11520 (4vr) 4 144 (4tt) 4 12 (4tt) 4 36 (4tt) 4 200 (4tt) 



a y 



The first line of this equation represents the one-loop result, the second line the two-loop 
result and the rest corresponds to the three-loop corrections. Since the Yukawa couplings 
enter the gauge j3 function starting from two-loop level only, it is enough to employ the 
one- loop RGEs for the Yukawa couplings for the precision we are aiming at. These can 
be derived in a similar manner from the general formula in ref. 
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The scalar self-couplings A, B and c that appear in eq. f|92|) only at the three-loop level 
are approximated as constants in our analysis by replacing them by their relations to the 
physical mass parameters M s , M 2 4, M Hc , Mx, and the gauge coupling a by using eqs. f )35|) 
and (I3"6"j) . The scalar self-coupling b that appears here can be approximated similarly by 
a constant using the SM Higgs mass SM and the mass of the W boson M w : 



3 2 ^H,SM 

4^ M 2 , 



(94) 
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